Solar Physics 
DOI: 10.1007/' 



Reflection and Ducting of Gravity Waves Inside the 
Sun 



K.B. MacGregor • T.M. Rogers 

© Springer •••• 

Abstract 

Internal-gravity waves excited by overshoot at the bottom of the convection 
zone can be influenced by rotation and by the strong toroidal magnetic field that 
is likely to be present in the solar tachocline. Using a simple Cartesian model, 
we show how waves with a vertical component of propagation can be reflected 
when traveling through a layer containing a horizontal magnetic field with a 
strength that varies with depth. This interaction can prevent a portion of the 
downward- traveling wave energy fiux from reaching the deep solar interior. If a 
highly refiecting magnetized layer is located some distance below the convection 
zone base, a duct or wave guide can be set up, wherein vertical propagation 
is restricted by successive refiections at the upper and lower boundaries. The 
presence of both upward- and downward-traveling disturbances inside the duct 
leads to the existence of a set of horizontally propagating modes that have 
significantly enhanced amplitudes. We point out that the helical structure of 
these waves makes them capable of generating an a-effect, and briefly consider 
the possibility that propagation in a shear of sufficient strength could lead to 
instability, the result of wave growth due to over-reflection. 

1. Introduction 

The ways in which internal-gravity waves can affect the compositional and dy- 
namical states of the solar radiative interior have received considerable attention 
in recent years. For example, wave-induced mixing of the layers just below the 
convection zone has been invoked to account for the observed depletion of lithium 
in solar/stellar photospheres (Garcia-Lopez and Spruit, 1991; Schatzman, 1996; 
Fritts, Vadas, and Andreassen, 1998), and it has been suggested that internal 
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waves can contribute to transport processes deep within the cores of the Sun and 
stars (Press, 1981; Press and Rybicki, 1981). The interaction between radiatively 
damped gravity waves and a mean shear flow, the forcing mechanism thought 
to be responsible for the observed quasi-biennial oscillation in the equatorial 
stratospheric layers of the Earth's atmosphere (see, e.g., Baldwin et ai, 2001, and 
references therein), has been investigated in the context of the time-dependent 
dynamics of the tachocline and the underlying stable region (Kumar, Talon, and 
Zahn, 1999; Kim and MacGregor, 2001; Rogers, MacGregor, and Glatzmaier, 
2008). There has also been much discussion concerning the effects of inwardly 
propagating waves on the overall internal solar rotation, particularly the long- 
term consequences of any angular momentum redistribution caused by waves 
for the rotational state of the radiative zone of the Sun (see Ringot, 1998, and 
references therein, for a summary). In this regard, it has been suggested (Talon, 
Kumar, and Zahn, 2002; Charbonnel and Talon, 2005; Talon and Charbonnel, 
2005; see also, Denissenkov, Pinsonneault, and MacGregor, 2008) that gravity 
wave interactions in the tachocline region and in the deeper interior can act in 
concert to produce near-uniform rotation of the core, in accord with helioseismic 
inferences (see, e.g., Charbonneau et ai, 1998). 

In a stably and continuously stratified fluid in which the effects of com- 
pressibility can be neglected, hydrodynamic gravity waves are generated when a 
localized region is perturbed by a small-amplitude disturbance that varies over 
a time scale that is longer than the period of adiabatic buoyancy oscillations 
at that position. On this basis, the inner and outer bounding surfaces of the 
Sun's convective envelope are likely sites of internal-wave emission into the 
contiguous, stable layers of the solar interior and atmosphere. These interfaces 
are deformed by convective fluid motions that overshoot and penetrate into 
adjacent stable regions where the local buoyancy period is typically short in 
comparison to the time scales characteristic of the perturbing flows. Kiraga et 
ai (2003), Rogers and Glatzmaier (2005a, 2005b) and Rogers, Glatzmaier and 
Jones (2006) have used detailed numerical simulations to study the excitation 
of gravity waves by overshooting plumes at the bottom of the Sun's convective 
envelope (r w 0.71 Rq). The inward propagation of the waves generated in this 
way must take them through the tachocline, which helioseismic analyses indicate 
has a central radius and thickness at equatorial latitudes of about 0.69 Rq and 
0.04 Rq, respectively (Charbonneau et ai, 1999). Numerous studies of flux tube 
formation and dynamics have delineated the links between the bipolar magnetic 
regions that are observed to emerge at the solar surface and magnetic flux stored 
in the subadiabatic layers beneath the convection zone (see, e.g., Schiissler, 1996; 
Fisher et al., 2000). Gollectively, the results of these investigations point toward 
the existence of a strong (« 10 — 100 kG), toroidal magnetic field within the 
tachocline, implying that the internal waves that traverse this region must be 
hydromagnetic rather than hydrodynamic in nature. 

The purpose of the present article is to examine one potential consequence of 
MHD modifications to the properties of gravity waves inside the Sun: namely, 
the possibility of reflection when a wave having a vertical component of propa- 
gation encounters magnetic conditions that change rapidly with depth below the 
convection zone. Some of the effects of a horizontal magnetic field at the bottom 
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of the convection zone on the propagation of internal gravity waves into the 
radiative interior have been investigated by Schatzman (1993). Since the length 
scale for radiative damping varies with wave frequency [uj] as {e-g-, Kim and 
MacGregor, 2001), consideration of the effects of reflection is likely to be most 
relevant for higher frequency waves that suffer less attenuation and are thus ca- 
pable, in principle, of propagating into the deep interior of the Sun (Press, 1981). 
Not only might a highly reflecting magnetic layer prevent some portion of the 
downward traveling energy flux from reaching the inner solar core, but it could 
contribute to the trapping of vertically propagating waves within a horizontal 
layer directly beneath the convection zone. That is, successive reflections of a 
wave from the base of the adiabatically stratified convection zone (wherein grav- 
ity waves are evanescent) and an underlying, depth- dependent toroidal magnetic 
field can lead to the formation of a duct or wave guide, effectively confining the 
wave to the region between the upper and lower reflecting surfaces, and allowing 
for horizontal propagation only. The existence and properties of such ducts for 
internal-gravity waves has been studied by atmospheric physicists in connection 
with the development and propagation of squall lines (Lindzen and Tung, 1976). 

In subsequent sections of this article, we develop preliminary descriptions of 
gravity wave reficction and ducting inside the Sun, using a simplified model 
for the magnetic and dynamical structure of the radiative layers comprised by 
the tachoclinc region. The present article is a companion to the recent articles 
by Rogers and MacGregor (2010, 2011), providing a conceptual basis for the 
interpretation of some of the detailed numerical results reported on in those 
articles. In Section 2, we derive the properties and propagation characteristics 
of MHD internal waves in a rotating, gravitationally stratified fiuid. The results 
of this section are then used in Section 3 to derive continuity conditions that 
enable us to treat the reflection of obliquely propagating gravity waves from 
a discontinuity in the strength of an otherwise uniform, horizontal magnetic 
field. Among other things, we derive an expression for the wave reflection coef- 
ficient and use it to deduce the combinations of wave frequencies and horizontal 
wavcnumbcrs for which the reficcting layer can function as the bottom of a 
duct. We obtain a solution for forced, hydromagnctic gravity waves in a duct 
of specified vertical extent, and examine the properties of the modes supported 
by such a structure when its lower bounding surface is perfectly reficcting. In 
the concluding section of the article (Section 4), wc summarize our results and 
briefly consider a few implications of gravity- wave reflection and ducting for the 
MHD physics of the solar interior, including: the helical structure of the waves 
and the possible relevance of this property to dynamo processes in the tachocline 
region, and the potential for over-refiection and instability in regions of strong 
shear. 

2. Gravity Waves in a Rotating, Magnetized Fluid 

When gravity is the only external force acting on an incompressible, stably 
stratified fiuid, propagating disturbances take the form of internal waves, driven 
by the fiuctuating pressure gradient and buoyancy forces (sec, e.g.^ Turner, 1973; 
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Lighthill, 1978). In the presence of rotation and a magnetic field, gravity- wave 
properties are modified through the additional influences of the Coriolis and 
Lorentz forces. To ascertain the nature of these changes, we consider infinites- 
imal, additive perturbations Sp, Sp, Su, and (SB to the density po, pressure po, 
velocity Uq, and magnetic field Bq of the mean equilibrium state. For an ideal 
MHD fiuid that rotates uniformly with a stationary angular velocity Jl, the 
behavior of these quantities is governed by the linearized equations 

V-(5u = 0, (1) 



^+^u.Vpo = 0, (2) 

Po —^ + 2po flxSu^ -VSp-gSp+—[ (V X Bo) x (SB + (V x 5B) x Bq ] , 
ot in 

(3) 

^ = V X ((Su X Bo) , (4) 

where g is the gravitational acceleration. The continuity Equation (1) expresses 
the assumed incompressibility of the fiuid, while Equation (2) indicates that 
fiuctuations in the density at any location are the result of the displacement of the 
mean density stratification by the fluctuating velocity field of the wave (Turner, 
1973; Lighthill, 1978). Equations (3) and (4) are, respectively, the linearized 
versions of the momentum equation (in the rotating frame of reference) and the 
induction equation for the inviscid, perfectly electrically conducting fluid. For 
simplicity, we have assumed Uq = 0, with the unperturbed equilibrium state 
determined by the balance between the pressure gradient, gravitational, and 
magnetic forces; the effect of a background flow flcld will be discussed in Section 
4. 

We arc ultimately interested in how the magnetic structure in the environs 
of the tachocline affects waves that propagate inward from the base of the 
convection zone. Since the radial extent of this region is estimated to be only 
w 10^^ Rq, we can facilitate the subsequent analysis by adopting a local Carte- 
sian coordinate system with x, y, and z-axes oriented parallel, respectively, to 
the directions of increasing 0, (^, and r in a spherical coordinate system with 
origin at the center of the Sun and polar axis aligned with the solar rotation 
axis. Relative to this Cartesian system, g = — g with g assumed constant, 
Bq = Bq(z) Gy, and n = U.^ Bx + with U,x = sin 9 and = 57 cos 9. 

The first two of these assumptions imply that density [po] and pressure [po] of 
the equilibrium state in the absence of waves are functions of z only. 

The treatment is further simplified by using two additional approximations: 
the Boussinesq approximation (Spiegel and Veronis, 1960; Turner, 1973), whereby| 
the variation in density is neglected in those terms in the momentum Equation 
(3) involving fluid inertia but retained in the buoyant force; and, the so-called /- 
plane approximation (e.(?.. Gill, 1982), under which the Coriolis force components 
arising from Vl^ in Equation (3) are dropped but those due to 0,^ are kept. The 
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first of these approximations requires that the vertical scale of the wave motion 
be smaller than the scale over which po varies, a condition which is well-satisfied 
in the radiative layers underlying the convection zone. The second approximation 
provides a reasonable description of wave dynamics for slowly rotating systems 
(i.e., when the rotation period is much longer than the period of adiabatic 
buoyancy oscillations), and when the domain of interest spans a narrow range in 
6 that is not located too near to the Equator. Although we confine our attention 
to disturbances that fulfill these constraints, it is also straightforward to carry 
out the derivation of wave properties presented herein without recourse to the 
/-plane approximation. 

Equations (1) - (4) admit of traveling wave solutions, having the general form 
5Q = 8Q{z) exp [i{kx + ly — ujt)], where SQ is any of the eight perturbation 
quantities Sp, Sp, Su, and (5B. Through a process of repeated substitution and 
elimination, each of these perturbations can be expressed in terms of the vertical- 
velocity fluctuation [Sur]. An equation for the amplitude [(5mz(z)] can then be 
obtained by substituting the resulting expressions for the perturbations into the 
z-componcnt of the momentum Equation (3), yielding 



_d 

d^ 



d 5u^ 



6u, 



0, (5) 



where ua — Bq/ ^AnpQ is the Alfvcn speed, / = 2 fi^ is the Coriolis parameter, 
and 



N = [ -(5/po) (dpo/dz) 



(6) 



is the Brunt-Vaisala or buoyancy frequency. 

As noted above, the remaining components of Equations (1)~(4) provide 
relations between 5uz and the other perturbation quantities. Specifically, the 
density fluctuation associated with the wave is given by 



6p^ Pq (iN^/guj) Sit;, 



(7) 



while the x and y-components of the velocity amplitude are 



5Ur 



iA:(a;2 - Pu\) - ojfl 
[k^ + P){uj^ - Pul) 



d Suz 
dz 



(8) 



and 



SUy = 



d Suz 
dz 



(k^ + P){uj^ - Pu\) 
The induction Equation (4) furnishes the means for determining S'B from Su 

I Sux 



(9) 



SB, = -Ba 



(10) 



SBy — Bo 



k Sux 



i_ d (Bg Suz) 
u) dz ' 



(11) 
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SB, 



-Bn 



I 5uz 



(12) 



with 5ux related to 5uz through Equation (8). Similarly, the total pressure 
perturbation [Jptot] is the sum of the fluid and magnetic pressure perturbations, 
and is expressible in terms of 5uz as 



Sptot = Sp + 



Bp 6B . 
47r 



d Si 



(13) 



For a medium in whieh both N and ua are independent of z, Equation (5) 
reduces to 



where 



d^ 5u, 
dz2 



+ Suz = 0, 



- 1 



1 - 



P - Pu\ 



n -1 



(14) 



(15) 



is the vertical component of the wavevector k = {k,l,m). In the ensuing dis- 
cussion, we focus on waves with horizontal components of propagation in the 
+ey-direction, so that k = (0, Z,m) with I > 0. Equation (15) can be recast in 
the form of an equation for uj, with solution 



P + TrP 



NH^ + fvPy Au\pPnP 



P + mP 



P + rrP 



1/2 



(16) 

in agreement with the dispersion relation given by Hide (1969), for the assumed 
directionalities of g, Bg, and k. Equation (16) contains as limiting cases a variety 
of wave modes, including Alfvcn waves (A^, / — > 0), hydrodynamic gravity waves 
(ma, / — > 0), inertial waves {N , ua — >■ 0), rotationally modified gravity waves 
{uA — >■ 0; e.g., Gill, 1982), MHD gravity waves (/ — >■ 0; e.g., Barnes, MacGregor, 
and Charbonncau, 1998), and MHD inertial waves (A'' — >■ 0; e.g., Lehnert, 1954). 

Insight into the propagation characteristics of MHD gravity waves in a ro- 
tating fluid can be gained by examining the behavior of the vertical wavevector 
component m as a function of lo. From Equation (15), it is apparent that for 
fixed values of the quantities N , ua, /, and I, m is real and vertical propagation 
is possible only for certain restricted ranges of uj values, as illustrated in Figure 
1. Specifically, inspection of the expression for mP indicates that mp > for 



U!f+<U!<U!B and u)f^ < uj < uja. 



(17) 



where 



UJB = {N^ + Pu\y/', ..;±=±i/+i(/2+4/\i)V2, coA = luA. (18) 

For all other values of uj, mP < and the waves are evanescent, decaying 
exponentially in the z-dircction. Figure 1 shows the regions of the ^-w-planc 
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Figure 1. The regions in the J-lj- plane for which gravity waves in a rotating, magnetized 
fluid can have a vertical component of propagation, as described in the text. The boundaries 
depicted correspond to rotation at a rate f /N = 10~^ (/ is the Coriolis parameter and N 
the buoyancy frequency) for waves with horizontal propagation in the +ey-dircction. When 
Bo, / = 0, vertically propagating hydrodynamic gravity waves are possible for any combination 
of I and u) for which u)/N < 1; that is, for the entire region below the horizontal dotted lino. 
For MHD gravity waves in the absence of rotation (/ = 0), vertical propagation is possible for 
waves within the region between the curve u) = ujg and the diagonal dotted line for u! = uij^. 
When the effects of rotation are included {/ ^ 0), vertical propagation is possible within the 
regions bounded by the curves uj = ujg, idf+ and ui = uij^, ijJf—- Rotation becomes increasingly 
important at low frequencies and long horizontal wavelengths, causing the curves for uj = uif-j- 
to diverge from each other. 

within which the inequahties (17) are satisfied for f /N = 10~^. This value is 
representative of conditions near the base of the convection zone at a latitude of 
about 30°, since 1] « 2.8 x 10"^ s^^ and iV w 2.5 x 10"^ s"\ the latter quantity 
estimated from the standard solar model of Bahcall and Pinsonneault (1995). 

By virtue of Equation (1) and the constraint V • (5B = 0, the waves under 
consideration are transverse, k ■ 6u ^ k ■ (5B = 0, so that both the velocity 
and magnetic fluctuations are contained in planes that are perpendicular to k. 
If /3 is the angle between a plane of constant phase and the vertical direction, 
then cos (3 = + m^)]^/^, where the variation of ni^/P with frequency 

can be deduced from Equation (15). Hence, near the upper boundary in Figure 
1 where a; ~ wb, the waves behave like ordinary hydrodynamic gravity waves 
for (Iua/N) <C 1, with m? /P sa 0, /3 w 0, and u) ~ N , implying horizontal 
phase propagation and nearly vertical fluid motions. For (Iua/N) 3> 1, w ss Iua, 
magnetic tension dominates the buoyant force, and the waves are like Alfven 
waves that propagate in the direction of the background magnetic flcld. As the 
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frequency is lowered at a fixed value of I the fluid motions become increasingly 
horizontal (/? « 7''/2), and buoyancy plays a smaller role in the wave dynamics. 
Near the line corresponding to a; = ujf^ in Figure 1, the waves satisfy the 
approximate dispersion relation 

with m? /P ^ 1. For {f /Iua) <SC 1, magnetic tension is again dominant and the 
waves arc Alfvcnic in character; when (f /Iua) > 1, the Coriolis force controls the 
fluid motions and the disturbances are much like hydrodynamic incrtial waves 
with u! f. At still lower frequencies, propagating waves (i.e., > 0) are 
possible within the region between the curves for to = ojA^^f-, wherein the 
approximate dispersion relation 

.^+[j^Y\.-fu\^Q. (20) 

holds. These waves have uj < Iua, and for {f /Iua) > 1 they are essentially the 
hydromagnetic incrtial waves described by Lehnert (1954) (see also Acheson and 
Hide, 1973). As uj uif-, the fluid motions take place in planes that approach 
a horizontal orientation, with the waves becoming evanescent for uj < ujf-. 

We note that like hydrodynamic gravity waves, the disturbances presently 
under consideration have vertical phase and group speeds, Vpz and Vgz, respec- 
tively, that are oppositely directed. In particular, by differentiating Equation 
(16), we obtain 

(21) 

from which it follows that for (say) m > 0, Vpz = Lom/n'^ > while Vgz < in 
both of the regions of vertical propagation delineated in Figure 1. 

3. Wave Reflection and Ducting 

We now use a simple model to investigate the reflection of vertically propagating 
MHD gravity waves in a layer containing a depth-dependent horizontal magnetic 
field. We consider small-amplitude, wave-like perturbations to a stationary, strat- 
ified medium in which the plane z = is a current sheet, separating the region 
(1) 2 > where Bq = Bi By from the region (2) z < where Bq — B2 By, 
with Bi and B2 constant fields. We follow the approach described by Fan (2001) 
to ensure that the density po of the unperturbed background atmosphere varies 
continuously with z, despite the prescribed jump in the strength of the magnetic 
field at z = 0. In particular, we express the background gas pressure and mass 
density in the form po = poa + Pqb, Po = Poo + Pob, where poo and poo are 
the pressure and density in an unmagnetized isothermal atmosphere, and pqb 
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and poB are the modifications to these quantities arising from the presence of 
the magnetic field Bq{z). Magneto- liydrostatic equilibrium of the background 
medium then requires that dpoo/dz = —gpoo and d[pQB + {BQ/8n)]/dz — —gpos- 
By choosing pos = 0, it follows that po is everywhere continuous and equal to 
the density of the unmagnetized atmosphere, po ~ poo- Integration of the second 
of the two equilibrium equations then yields pos = {Bf — i?|)/87r, from which 
it can be seen that although a jump in the gas pressure exists at z = 0, the 
total pressure, ptot = Po + Bq/Stt, is continuous there. From the ideal-gas law, 
the discontinuity in po is associated with a jump Tqb ~ Too [pob /poo) in the gas 
temperature at the current sheet location. The fractional changes in p and T are 
both w /3^^ and are small since f3 = {Sttpoq/Bq) « 10^ for physical conditions 
like those at the base of the convection zone. Note that because the effects of 
compressibility have not been included in the analysis of Section 2, the sound 
speed (i.e., Tq) does not appear in the dispersion relation (15), making that 
result applicable to the description of waves in either of the regions 1 and 2 
deUneated above. 

We emphasize that the adopted background atmosphere model is intended as 
an approximate representation of a neutrally buoyant, magnetic layer in the ra- 
diative region beneath the convection zone. We omit (among other things) treat- 
ment of the the energy balance within the layer, assuming instead a piecewise- 
isothermal temperature distribution in the unperturbed medium that is consis- 
tent with the force balance that prevails therein. This simplified model does, 
however, allow us to isolate and explore magnetic effects on internal wave prop- 
agation under solar-like conditions since, for po continuous and A'' assumed 
constant. Equation (15) indicates that a wave traveling from region 1 to region 2 
will be affected solely by the discontinuity in Bq at z = 0. Acheson (1976) utilized 
a similar equilibrium model to investigate the over-reflection of hydromagnetic 
gravity waves in a medium that also contains a strong shear flow. An alternative 
treatment of the discontinuity in po implied by the assumed Bo{z), namely, a 
jump in Po with To continuous at z = is discussed later in this section; in this 
case, wave propagation is affected by both the variations in Bq and po. 

We focus on a plane wave of frequency w that originates in region 1 with 
K = (0, /, m), where / and m are both > 0. Such a wave has its vertical component 
of propagation (i.e., Vgz) in the direction of decreasing z (see Equation [21], and 
the discussion thereof); in traveling downward from region 1 into region 2, the 
wave will encounter a discontinuous change in the Alfven speed, from uai to 
UA2, at z = 0. We anticipate that, in general, both incident and reflected waves 
will be present in region 1, while region 2 will contain a transmitted wave. The 
vertical velocity fluctuation arising from each of these components is a plane 
wave of the form Su^ = Su^^ exp [i{ly — ujt)], where diiz satisfies Equation (5). For 
the conditions outlined in the preceding paragraph, it is clear that the relevant 
solutions to Equation (5) are 



Suz{z) = 6uj exp (imiz) -I- Sur exp (— imiz). 



(22) 



in region 1, and 



Suziz) = 5ut exp (iTO2z), 



(23) 
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in region 2, where nii (TO2) is the vertical wavevector component (see eq. [15]) 
in the upper (lower) region, and diii, 5ur, and 5ut are the constant amplitudes 
of the incident, reflected, and transmitted waves at z = 0. 

The wave solutions in the upper and lower halves of the domain can be 
connected across the interface through the application of physical conditions 
expressing the continuity of the vertical velocity and total pressure perturbations 
at z = 0, 



5uz cxp [i{ly — Lijt)] 



0, 



Sptot exp [i{ly - Lot) 



0. 



(24) 



For Sptot given by Equation (13) with po assumed continuous at the interface, 
it follows directly from the imposition of these constraints that u! and I are the 
same for the incident, reflected, and transmitted waves, and that 



where 



6ui 



m2 
Q= 

nil 



Sill, Sut 



Sill 



'■A2 



/2,,2 



l-[uj,f/{u:^-Pu\^) 
1 - [ ufliLO^ - Pu\^) 



(25) 



(26) 



A quantity of particular relevance to the present investigation is the reflection 
coefficient, i?, defined as 



R = 



5ur 




1-9 


Suj 




i + g 



(27) 



using the first of Equations (25). 

A qualitative understanding of the reflection of MHD gravity waves can be 
developed by considering the ranges of I and uj for which vertical propagation is 
possible in regions 1 and 2. The curves in Figure 2 delimit the portions of the luj- 
plane in which mi (solid lines) and m2 (dashed lines) are real, for the particular 
parameter values UA2/UA1 = 5.0, 0.2 and f /N = 10~^. In the following discus- 
sion, we focus on the particular case depicted in panel (a) for [Iuai/N) — 10^^; 
analogous considerations can be applied to determine the propagation character- 
istics of waves corresponding to other values of / and ua2/uai - In proceeding from 
high to low frequencies along the vertical dotted line in the figure, the discussion 
of Section 2 indicates that propagating waves can exist in both regions 1 and 2 
when ujfj^2 < uj < lobi ~ ujb2, where the frequencies defining this interval are 
given in Equations (18). In the range uja2 < uj < w/+2, waves can propagate in 
region 1 but are evanescent in region 2, while for w/-2 < uj < uja2, propagation 
is again possible in both regions. For uj < w/-2, wave propagation is precluded 
in region 2, while waves in region 1 are traveling for ujf^i < uj < ujf-2 and 
ujf-i < UJ < uJAi but evanescent for ujai < uj < i^/+i and uj < ujf^i. 

The information obtained from Figure 2 concerning wave propagation on 
either side of the current sheet at z = can be used to infer something about the 
behavior of the reflection coeffcient i? as a function of I and uj. For waves that can 
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Figure 2. MHD gravity-wave propagation in a medium with Alfven speed Uj\\ for z > 
and uj^2 for ^ < 0. The curves are analogous to the boundaries identifled in Figure 1, and 
deUneato the portions of the /-Lj-plane whore waves have m? > in the upper (solid lines) or 
lower (dashed lines) half of the domain, as explained in the text. The results shown in the two 
panels were obtained for the parameter values f /N = 10^^ with ua2/uai = 5.0 (panel a) and 
0.2 (panel b], 

propagate in region 1 but are evanescent in region 2, < 0, and Equations (26) 
and (27) indicate that q is imaginary and R = 1. In this case, the discontinuity in 
UA at the interface between regions 1 and 2 acts like a perfect reflector. For values 
of I and uj such that vertical propagation is possible in both halves of the domain, 
q is real, R < 1, and the jump in ua is partially reflecting. These deductions 
are verified by inspection of Figure 3, wherein we show the variation of R with 
w for UA2/UA1 =5.0 and 0.2, and {Iuai/N) = 5 x lO^^, lO^^, lO^^, 5 x lO^^; 
these values of the horizontal wavevector component are indicated by the vertical 
dotted lines in Figure 2. For reference, note that if the mass density at the base of 
the convection zone is taken to be about 0.2 g cm~^, then uai ~ lO'* cm s^^ for 
a field strength Bi in the range 10 - 20 kG. With N w lO'^, (Iuai/N) = 10"^ 
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Figure 3. The reflection of MHD gravity waves in tiie composite medium of Figure 2. The left- 
hand (right-hand) panels pertain to wave propagation in a medium with ua2I'>'-A\ = 5.0 (0.2). 
The reflection coefficient R, defined as the ratio of the vertical velocity amplitude of the 
reflected wave at 2: = to that of the wave incident there, is shown as a function of frequency 
for waves with {Iuai/N) = 5 x 10~* (panels a and b), 10"'^ (panels e and d), 10"'^ (panels 
e and f), and 5 X lO"'^ (panels g and h); these values correspond to the vertical dotted lines 
in panels (a) and (b) of Figure 2. The long vertical ticks at the top of each panel indicate the 
frequencies i^j—, (^Ay ^f+y s-^d uib (sec Equations (17) and (18) and the discussion thereof) 
that demarcate the frequency intervals within which wave propagation is possible for the given 
I value in region 1 (ua = UAi) or region 2 {ua = UA2)- 
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then implies that I « 10^^ cm^^ « 2tt/Hp, where Hp 0.08 Rq) is the 
pressure scale height at the convection zone bottom. 

Referring to the discussion of the preceding paragraph and to Figure 2, it is 
readily seen from Figure 3 that within those frequency intervals for which ml > 
and 7712 < 0, i? = 1. Likewise, at frequencies for which freely propagating waves 
are possible in both regions, R < 1, with R undefined at frequencies for which the 
waves are everywhere evanescent. As is apparent from the figure, as (Iuai/N) 
increases, the frequency intervals with i? = 1 shift to higher uj/N, while the 
intervals with R < 1 narrow. For downward propagation in a medium in which 
the strength of the horizontal field increases with depth, there are frequency 
intervals having i? = 1 for all values of (Iuai/N); in the case of a horizontal 
field strength that decreases with depth, perfect reflection occurs only for low 
frequency MHD inertial waves having {Iuai/N) <^ 1. For frequencies at which 
the magnetic layer is partially refiecting, i? is a strongly varying function of uj, 
with the interface becoming perfectly transmitting (i.e., (7=1 and i? = 0) for 
waves with {co/N)"^ « ^{1 + {Iuai/N)^ [1 + {ua2/uai)^ ]}■ Behavior analogous to 
the cases shown in Figure 3 is seen for other values of the ratio ua2 / UAi , the only 
differences being in the extents of the frequency ranges for which the reflection is 
perfect or partial. For the special case in which uai ~ 0, ua2 0, propagation in 
region 1 is possible for { = j) < to < ujb { — N). Hence, for {Iuai/N) large 
enough that ujf-2 > /, the variation of R with uj for uj > is as depicted 

for the cases with ua2/uai > 1 in Figure 3, with R = 1 for w/-2 > > /• 
For [luAi/N) such that a;/_2 < /, the lowest frequency propagating waves in 
region 1 are partially reflected if w/-2 < / < '^A2 and perfectly reflected if 

^A2 < f < W/+2. 

The density and total pressure in the background equilibrium atmosphere 
considered throughout this section are continuous across the current sheet at 
which the magnetic field, gas pressure, and temperature change discontinuously. 
For completeness, we note that the jump in gas pressure required to maintain 
equilibrium in the presence of the jump in Bq could likewise have been provided 
by a jump in density with the temperature and total pressure continuous. For 
wave reflection in this case, with the density change across the sheet Ap = p2 — Pi 
assumed small {Ap/po <^ 1), the linearized continuity condition corresponding 
to the second of Equations (24) becomes (see, e.g., McKenzie, 1972; Delisi and 
Orlanski, 1975) [ Sptot — 9Po£, ]i = 0, where ^ = iSiiz/uj is the amplitude of the 
wave-induced vertical displacement of the z = surface. Following a procedure 
analogous to that used in deriving the reflection coefficient (27), we find 



R 



l + q + 

where q is still given by Equation (26) and 



(28) 



^ %o mi {Lu-^-Pu\^)^'^-uj^p- ^^^^ 

Note that for frequencies corresponding to waves that can propagate in region 
1, TOi and Q cx {Ap/ pq) cx I3~^ ^ 1 are both real quantities. For frequencies 
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corresponding to waves that are evanescent in region 2, ?7i2 and q are imaginary, 
and inspection of Equation (29) reveals tliat i? = 1 as was the case for reflection 
from a current sheet at which the density was continuous. For frequencies corre- 
sponding to waves that can propagate in both regions 1 and 2, g is real and > 0, 
and it is readily shown that i? < 1; as before, such waves are partially reflected. 
More specifically, in this case we find 

'-'^%ii'= i;-'':+g: <i, ,30) 



l + qj (l + g)2 + Q2 

indicating that for frequencies at which partial refiection of the incident wave 
occurs, the reflection coefficient coefficient is enhanced by a small amount relative 
to the R value that obtains for Ap — (c/. Equation (27)). 

As described in Section 1, gravity waves are excited at and travel downward 
from the bottom of the convection zone, any upward propagation prohibited by 
the super-adiabatic stratification that prevails in the region above. Hence, re- 
flection from a somewhat deeper lying magnetized layer could effectively confine 
the vertical propagation of waves having i? w 1 to a thin slab-like region below 
the convection zone. To investigate the circumstances under which magnetic 
structure might contribute to the formation of such a duct or wave guide in the 
outermost layers of the radiative interior, we consider a configuration in which 
the strength of a y-directed magnetic field changes from Bi to B2 at depth 
z ~ —d beneath the convection zone base (z = 0). Specifically, we identify the 
layers —d < z < and z < —d with the regions 1 and 2, respectively, of the 
previous discussion, and employ the solutions given in Equations (22) and (23) 
to describe the waves present in the upper and lower portions of the domain. 
As before, the continuity of the solutions at 2 = —d is ensured by applying the 
conditions (24), and we utilize the requirement that Suz{Q) ~ Suq to account for 
the excitation of waves at z = 0. 

Following the procedure outlined above, we find that the vertical component 
of the wave velocity is 

exp [imi(z -I- d)] -I- Cfl exp [-imi(z -I- d)] 

du^ = duo T. T— -7^ r~- T\ ' (31) 

exp (imid) + Cr exp (— imid) 



in region 1, and 



c Ct exp [im2{z + d)] 

ou, = duo — , (32) 

exp (imid) + Cr exp (— imid) 



in region 2, where the coefficients Cr and Ct are 

Cfl = "TT^' C't = — — , (33) 
1 + q 1 + g 

with q still given by Equation (26). The effectiveness of the layer — d < z < 
(i.e., region 1) in functioning as a duct for gravity waves depends upon the 
magnitude of the refiection coefficient [R] of the magnetic layer that serves as 
its lower boundary. For specificity, in the following we consider the case of wave 
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propagation into a region of increasing horizontal field strength, as depicted in 
Figure 2a. On the basis of the results depicted in Figure 3, we can expect to 
trap little of the wave energy in the frequency band ujf-2 < uj < ujbi, since 
R = \ only in the narrow interval wa2 < w < w/+2 with R everywhere else. 
The reflectivity of the bottom layer is considerably higher at lower frequencies, 
however, since all waves with uj < w/-2 that can propagate in region 1 are 
evanescent in region 2 and have R= 1. We therefore expect that the duct that 
is formed when region 1 is sandwiched between the convection-zone base on the 
upper side and region 2 below will suffer minimal leakage for propagating waves 
with < oj < w/_2. 

These expectations are confirmed through examination of the amplitude of the 
total pressure perturbation Sptot associated with MHD gravity waves in region 
1. Evaluating (5j5tot using Equation (13) with the solution given in (31), we find 



5ptot = - 



po Suq mi 



1 - 



'Al 



{Pr + iP^) , (34) 



where 



Pr 



D 



(1 — R^) cos (miz). 



(35) 



P 



D 



{(1 + R^) sin (miz) + 2C]^ sin [mi{z - 



-2d)] -2CJj cos [mi(z + 2( 



and 



D = l + R^ + 2Q cos (2mid) + 2Q sin {2rnid). 



(36) 



(37) 



In Equations (34) -(37), the reflection coefficient is defined as i? = | Cr |, with 
Cr = C}j + i C]^. In Figure 4, we show the frequency dependence of the pressure 
fluctuation | 5ptot \ at the tops (i.e., z = 0) of ducts with vertical thicknesses 
Id = 10~^ and Id — \. For the purpose of comparison with previous figures, 
the present results pertain to waves with (Iuai/N) = 10~^ in a medium with 
//TV = 10~^ and ua2/uai — 5.0. For the parameter values assumed throughout 
this discussion, the implied thickness of the ducting region when Id ~ l\s d ~ 10^ 
cm w 10-2 

As is apparent from the figure, at higher frequencies where i? ^ 1, the waves 
are little affected by the change in magnetic conditions at the interface between 
regions 1 and 2 and propagate nearly freely. From the discussion of Section 
2, for the adopted parameter values, these waves behave like hydrodynamic 
gravity waves, with | Sptot \~ Po Suq ujuii/l'^ « po Suq N/l. However, within the 
frequency intervals LOA2 < < w/+2, < uj < ojf-2, and ujf-i < uj < cuai, 

R = 1, and the duct supports a sequence of wave modes, identifiable in the figure 
as significant enhancements in the value of | Sptot \ at discrete frequencies. The 
restoring force for these modes is a combination of the pressure gradient, Coriolis 
and magnetic forces, the exact balance depending upon the magnitude of the 
quantity (f/luAi) (see Equations (19) and (20), and the discussion thereof). 
The vertical propagation of the waves is restricted to the layer —d < z < 
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Figure 4. The amplitude of the total pressure perturbation (solid line; see Equation (13])) 
at 2 = for MHD gravity waves in a duct with vertical extent d, as discussed in Section 
3 of the text. The results depicted were obtained for {Iuai/N) = 10"^, f/N = IQ-^, and 
*'A2/*'Ai = 5.0, with Id = 10~^ in panel a and W = 1 in panel b. The dashed line in both 
panels represents the reflection cocffcicnt [_R], while the dash-dotted line in panel (a) shows 
the value of the denominator [D] given by Equation (37) . The vertical dotted lines mark the 
frequencies given by Equations (18), evaluated in regions 1 and 2 (see also Figure 3). 
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by virtue of the structure of the duct, but they are able to propagate without 
attenuation in the horizontal direction. A similar behavior is seen in models 
for ducted hydrodynamic gravity waves in the Earth's atmosphere (see, e.g., 
Lindzen and Tung, 1976). 

In the case under consideration, the peaks at which the total pressure per- 
turbation is intensified are produced when waves in the layer that have opposite 
senses of vertical propagation constructively interfere. Inspection of panel a in 
Figure 4 reveals that the peaks coincide with the zeros of the denominator 
defined in Equation (37); these occur when 

2mid - 6* = (2?i + 1) TT, n = 0,1,2,... (38) 

where 9 — tan^^ (C^r/^b)- Since rrii increases with decreasing frequency (see 
Equation (15)), this relation yields a set of mode frequencies a;„ that decrease 
as n gets larger. At fixed I, the horizontal phase speeds of these wave modes, 
Vpy = (w„/Z)/[l + (m/l)^], become slower for increasing n. Condition (35) also 
establishes the vertical structure of the modes in the duct; for {6 /Air) ^ 1, the 
wave with frequency w„ has (2n-|- 1)/4 vertical wavelengths within the duct width 
d. Taken together, these properties suggest that low-frequency, high-n modes are 
likely susceptible to radiative, viscous, or resistive dissipation. The peaks seen 
in Figure 4 occur for frequencies low enough that i? = 1; this means that the 
peak with the highest frequency has n = 3 in the case Id = 10^^ and n = 13 in 
the case /d = 1. In this latter example, the peaks corresponding to lower values 
of n are evident in the figure, albeit with reduced amplitudes resulting from the 
fact that the partially refiecting (i.e., R < 1) lower boundary allows a fraction 
of the vertically propagating wave fiux to escape from the duct. 

It was pointed out in Section 2 that MHD internal waves are transverse, the 
fluid motions associated with them taking place in planes perpendicular to the 
direction of phase propagation. Due to the influence of the Coriolis force, Su 
does not maintain a flxed orientation within a given plane but instead rotates, 
giving the wave a helical structure (see, e.g., Moffatt, 1978). Waves of this kind 
may have some relevance to dynamo processes inside the Sun, since helical fluid 
motions are required in order to produce an a-effect. The kinetic helicity of a 
single wave can be evaluated by computing the scalar product of the velocity fluc- 
tuation Su with the vorticity V x 6u. Using Equations (8) and (9) to determine 
Su for a freely propagating wave with Suz = Suq exp (iV'), ip = (ly + rniz — cot), 
we flnd 

Suq — ^ cos ip, Suz — Suq cos ij}, 

(39) 

nK = Su-{y xSu) = -miSul {-^^^L_^ + (40) 

for the kinetic helicity. Similarly, for MHD gravity waves in a magnetically 
defined duct, we use Equation (31) for Siiz to obtain (after considerable ma- 



Suq 



oof 



/2,,2 



— ^ sm %}, Suy 



yielding 
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Figure 5. The absolute magnitude of tlie liinetic helicity of MHD gravity waves for 
{luAi/N) = 10~^ and f/N = lO""^. The dashed and solid lines represent the helicities 
associated with a wave that propagates without reflection (Equation (40)) and a ducted wave 
for the case Id = 1 (Equation (41)), respectively. In the latter case, the lower boundary of the 
duct was taken to be layer with «A2/f Ai = 5.0. The insert provides additional detail regarding 
the behavior of T-Lk in the vicinity of the frequencies ti^/+2i '^A2i '^/-2- 



nipulation) 

where R is the reflection coeffcient and D is given by Equation (37). 

In Figure 5, we show the absolute magnitudes of the kinetic hehcities asso- 
ciated with both freely propagating and ducted waves, for {Iuai/N) = 10^'^ 
and f/N = 10""^, and for a duct with thickness Id = 1; the results for Id — 
10^^ are quite similar to those shown in the figure. From Equations (40) and 
(41). it is clear that when / (= 2n cos 9) > 0, the helicity of waves with 
TOi > that transport energy downward from the upper boundary is nega- 
tive for Lj > luAi and positive for cj < Iuai- The former of these frequency 
domains includes waves for which buoyancy contributes to the wave motion 
(a;/+i < w < wsi), while the latter corresponds to the hydromagnetic in- 
ertial waves (w/-i < < ujai)- For decreasing w in both frequency bands, 
m/l becomes 3> 1, leading to wave- induced fluid motions that are increasingly 
horizontal and growing helicity. An example of this behavior is provided by 
gravity waves in the frequency interval 5 x 10^^ < oj/N < 5 x 10^^, for which 
the 'Hk/{15uq) sa —{f/bj){mi/l)^ « —fN^/co'^; at lower frequencies, for the 
hydromagnetic inertial waves, 'Hk/{15uq) ~ {i^/^Y^ where n = {l"^ + m^)^/^ is 
the total wavevector. 
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As is evident in the figure, the presence or absence of a reflecting interface 
has little impact on the wave helicity for uJf+2 < io < lobi where R < 1. 
Apart from the enhanced helicity of ducted waves with frequencies near the 
lower end of this interval, the resuhs for both cases depicted in the figure are 
nearly indistinguishable from one another. At lower frequencies, however, where 
waves are evanescent in region 2 and efficient refiection leads to the existence 
of the modes pictured in Figure 4, the hehcity of the ducted waves vanishes. 
According to Equation (40), for waves with < w < w/-2, downward 

(mi > 0) propagating waves have negative helicity and upward propagating 
waves (mi < 0) have positive hehcity. Hence, when R = 1, the presence of equal 
fluxes of waves with opposite helicities in the duct yields zero net helicity, as 
indicated by Equation (41). 

Because of their helical structure, it is of interest to ascertain whether these 
waves can contribute to the generation of a magnetic field through the production 
of an a-effcct. To do this, we follow Moffatt (1978) and consider the case in which 
wave propagation takes place within a weakly dissipativc medium having mag- 
netic diffusivity 77. Using the fiuctuating magnetic and velocity fields associated 
with a single freely propagating wave to calculate the mean electromotive force 
£ = {Su X 6'B) = q;Bo, we find 



As indicated by inspection of Figure 5, this should also be a reasonable esti- 
mate of the a-effect arising from waves that undergo partial reflection from the 
magnetic layer in region 2. Adopting rj = 10^ cm^ for the diffusivity at 
the convection zone base and assuming Hk / (ISuq) w 10^ — 10^ (Figure 5) with 
I « 10~^ cm~^ and Suq ~ few x 10^ cm s~^ for uj/N w 5 x 10~^, we obtain 



4. Summary and Discussion 

We have elucidated the physical properties of the internal waves that are likely 
to be present in the stable layers underlying the solar convection zone. We have 
used a simple, Cartesian model to conduct an exploratory examination of the 
reflection of vertically propagating waves in a region containing a horizontal (i.e., 
toroidal) magnetic field whose strength varies with depth; further investigation 
of these effects will require the use of a more realistic computational model (see, 
6.17., Rogers and MacGregor, 2010, 2011). For given values of I and uj, waves that 
travel downward from the base of the convection zone can undergo reflection with 
i? = 1 if they encounter a layer in which, because of the changing background 
magnetic conditions, they become evanescent. The presence of an effectively 
reflecting magnetized layer below the wave source region at the bottom of the 
convection zone can lead to the formation of a duct or wave guide, a structure 
that limits the vertical propagation of the perturbations and supports a set of 
horizontally propagating modes with enhanced amplitudes. In this regard, we 




(42) 
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note that inspection of the results for R depicted in Figure 3 indicates that 
for each value of the horizontal wavevector, there is a corresponding range of 
wave frequencies for which i? « 1. If the energy flux of gravity waves excited 
by convectivc overshoot has the relatively flat distribution in frequency seen in 
the simulations of Kiraga et al. (2003) and Rogers and Glatzmaier (2005), then 
reflection could prevent a significant fraction of the emitted wave energy from 
reaching the deep solar interior. Such redirection and trapping of inward traveling 
waves would no doubt have implications for models in which the helioseismically 
inferred near- uniform rotation of the radiative interior is a long-term consequence 
of angular momentum redistribution by gravity waves {e.g., Charbonnel and 
Talon, 2005, and references therein). 

The internal waves studied herein have helical structure, raising the possibility 
that such disturbances could contribute to the amplification of magnetic fields 
through dynamo action. Along these lines, Schmitt (1984, 1987) has investigated 
a model in which unstable magnetostrophic waves, driven by magnetic buoyancy, 
provide the a-cffect for a dynamo located at the bottom of the convection zone. 
In the case considered herein, although the estimated a-effect produced by a 
single freely propagating wave is small, the cumulative effect of a superposition of 
waves may be larger. The results of Figure 5 suggest that a further enhancement 
of the net helicity (and thus, of a) of waves with uj w might also be achieved 
through propagation in a leaky duct with a partially reflecting bottom. Ducted 
modes for R = 1, however, exhibit zero net helicity, a result of the cancellation of 
the positive and negative helicities of upward and downward traveling waves. In 
addition, mixing processes driven by these modes could affect the abundances of 
light elements in the outer layers of the radiative interior, although a quantitative 
assessment of these and related effects requires a more realistic representation 
of the structure of this region, together with treatment of such influences as the 
radiative damping of low frequency waves and irregularly shaped, non-horizontal 
reflecting surfaces {e.g., Phillips, 1963). A duct thickness d w 10^ cm below 
the convection-zone base represents an appreciable fraction of the distance over 
which mixing must occur in order to ensure destruction of lithium by nuclear 
processes. 

Finally, the analysis presented in the preceding sections did not directly ac- 
count for the rotational shear flow that is the salient dynamical feature of the 
tachocline region. Note that in the presence of a mean background flow of the 
form u = u{z) By, the vertical components of the wavevector and group velocity 
become (see, e.g., Barnes, MacGregor, and Charbonneau, 1998) 



and 



(cj - luY - Pu\ 



~ ^{uj-Iu)k^' 



1 



(43) 



(44) 



respectively, where k = (0, 1, m) and rotation has been neglected {i.e., / — 0), for 
simplicity. We consider the reflection of a wave that is incident on the interface 
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between the region z > where ua = uai and u = ui, and the region z < 
where ua — ua2 and u ~ U2, adopting the configuration and nomenclature of 
Section 3. Utilizing the procedure given in Section 3 to derive the reflection 
coefficient, we again obtain R=\ (1 — g)/(l + q) |, but with 



instead of Equation (26). 

For conditions such that waves are propagating in region 1 but evanescent in 
region 2, we anticipate that wave reflection should, apart from Doppler shifts 
arising from the advection of waves by the flow, qualitatively resemble the re- 
sults obtained in Section 3 assuming u = 0. In particular, such a magnetized 
shear layer should support horizontally propagating, ducted modes of the kind 
investigated in Section 3, when the effect of an overlying convective region on 
the vertical propagation of waves is accounted for. Alternatively, note that if. 



then by Equation (43) m^, rri^ > and vertical propagation is possible in both 
regions 1 and 2. For the case in which ui < lj/I < U2, Equation (44) indicates 
that Vgz < for the incident wave (mi > 0) in region 1 but Vgz > for the 
transmitted wave in region 2, unless 7712 is chosen to be < 0. With this choice, 
however, it follows from (46) that g < and R> 1, so that the amplitude of the 
reflected wave exceeds that of the wave incident on the interface. In this case, the 
wave has undergone "over-reflection" , with the amplitude of the reflected distur- 
bance increased through interaction with the specified background shear flow (see 
Acheson 1976, and references therein). From the preceding analysis, we conclude 
that the occurrence of this process requires i) the existence of a strong shear (in 
fact, the Richardson number Ri = [N / {du / dz)]'^ must be Y4i see Acheson, 
1976), ii) that the horizontal flow speed [u] somewhere in the shear layer exceed 
uj/l, and in) that vertical propagation be possible on both sides of the interface 
{nil, "^2 > 0)- Since the vertical shear within the tachocline region is thought to 
be characterized by Ri ^ Y4 (Schatzman, Zahn, and Morel, 2000), it is unlikely 
that over-reflection of MHD gravity waves takes place. However, if the prevailing 
shear properties were to be conducive to the occurrence of over-reflection, then 
the magnetic shear layer would necessarily become unstable since, for waves in 
the appropriate frequency range, each successive reflection from the interface 
would increase the wave amplitude by a factor of R (> 1; see also Acheson, 
1976). Among the consequences of over-reflection of internal waves from the 
lower boundary of such a ducting region would be a non-zero net helicity for the 
modes in the presence of rotation (i.e., / 7^ 0), since the differing amplitudes of 
the waves traveling in the ie^-directions implies that their respective helicities 
will be unequal in magnitude. Under these circumstances, it might be possible 
for the layer to function as a dynamo, amplifying the field through fluid motions 
whose energy source is the ambient shear flow. 



q 



1712 \ (g; - lu2)'^ - Pu\^ 



(45) 



mi [uj — lui)^ — Pu^ 



(46) 
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